Influence of Fermi surface topology 
on the quasiparticle spectrum in the vortex state 
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We study the influence of Fermi surface topology on the quasiparticle density of states in the 
vortex state of type II superconductors. We observe that the field dependence and the shape of 
the momentum and spatially averaged density of states is affected significantly by the topology of 
the Fermi surface. We show that this behavior can be understood in terms of characteristic Fermi 
surface functions and that an important role is played by the number of points on the Fermi surface 
at which the Fermi velocity is directed parallel to the magnetic field. A critical comparison is made 
with a broadened BCS type density of states, that has been used frequently in analysis of tunneling 
data. We suggest a new formula as a replacement for the broadened BCS model for the special 
case of a cylindrical Fermi surface. We apply our results to the two gap superconductor MgB2 and 
show that in this particular case the field dependence of the partial densities of states of the two 
gaps behaves very differently due to the different topologies of the corresponding Fermi surfaces, in 
qualitative agreement with recent tunneling experiments. 

PACS numbers: 74.25.Op, 71.18.+y , 74.50.+r, 74.70.Ad 



I. INTRODUCTION 

The quasiparticle density of states is one of the char- 
acteristic properties of a superconductor. It contains in- 
formation about the gap function and can be probed by 
various experimental techniques such as tunneling spec- 
troscopy, photocmission, specific heat, nuclear magnetic 
resonance (NMR) etc. These techniques are often used 
to extract information about the superconducting gap in 
a given system. Recently, there is growing interest to 
study the quasiparticle excitations in the vortex state of 
type II superconductors, particularly in unconventional 
superconductors ii»2iMiMii For example, the zero energy 
density of states in a d-wave superconductor is expected 
to vary as the square-root of the magnetic field B instead 
of the linear B-field variation expected in a dirty s-wave 
superconductoriSi2ii£ However, the relation between the 
gap function and the quasiparticle density of states is 
less direct in the vortex state than in the Meissner state. 
The reason is that there are important modifications of 
the density of states due to bound states in the vortex 
cores and the supercurrents running around the vortices, 
particularly at higher magnetic fields. In this case quasi- 
particles are excited with respect to the local supercur- 
rent flow and their energy appears 'Doppler-shifted'pii 
In analysis of tunneling spectra these effects are often 
modeled as a 'smearing effect' of the zero field density 
of states4^*i However, as we will show below, the quasi- 
particle excitations above the moving supercurrent con- 
tain information about the momentum dependent Fermi 
velocity, especially its direction relative to the applied 
magnetic field. For this reason the particular shape of 
the Fermi surface has an important influence on the mo- 
mentum averaged quasiparticle density of states in the 
vortex state. 

Here, we want to present a detailed study of this in- 





FIG. 1: Different shapes of the Fermi surface that are investi- 
gated within this work. In the upper row we have on the left 
hand side a spherical Fermi surface and on the right hand side 
an elliptical Fermi surface. In the lower row we show on the 
left hand side a distorted cylinder with a distortion parameter 
e c = 0.163 and on the right hand side a half-torus. 



fluence of the Fermi surface structure on the quasiparti- 
cle density of states in the vortex state. We are going 
to compare the different Fermi surface types shown in 
Fig. ^ spherical and elliptical Fermi surfaces as they 
are used commonly in isotropic and anisotropic effective 
mass models, a cylindrical Fermi surface with a small c- 
axis dispersion as appropriate for layered systems, and a 
Fermi surface of half-torus shape which is relevant for the 
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7r-band in the recently discussed superconductor MgB2, 
for example^ Our calculations are based on quasiclassi- 
cal Eilenberger theory using a method that was originally 
introduced by Peseta to describe type II superconduc- 
tors at high magnetic fields and that we have generalized 
to unconventional pairing symmetry recently^ In that 
work we have shown that this method leads to much more 
accurate results than the frequently used 'Doppler-shift 
method', particularly at higher magnetic fields, mainly 
because it takes into account the contributions from vor- 
tex core states properly. Our results presented below 
show that the structure of the Fermi surface affects the 
density of states in the presence of a magnetic field quite 
dramatically. We are going to describe our calculational 
procedure in section [HI In section PHI we introduce the 
different Fermi surface types used in the present study. 
Section IIVI contains our results and interpretation in 
terms of characteristic Fermi surface functions. In sec- 
tion[V]we compare our results with the model of a broad- 
ened BCS type density of states that has been used fre- 
quently for analysis of experimental data. In section IVTI 
we study the particularly interesting case of the two gap 
superconductor MgP>2, which possesses four Fermi sur- 
faces of two different types: two cylindrical c-bands and 
two 7r-bands, which can be modeled by the half-torus 
in Fig. E 12 These different Fermi surface types lead to 
very different field dependencies of the partial densities 
of states for the two gaps. 
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where x and y are the coordinates in the plane perpendic- 
ular to the magnetic field B and the complex quantities 
loi and u>2 span the vortex lattice A. J\f is a normal- 
ization factor that has to be chosen in such a way that 
\ip\\ 2 averaged over a unit cell of A equals unity. The 
normalized local density of states can be obtained from 
the real part of the analytical continuation of g, averaged 
over the Fermi surface 
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denotes an average over the Fermi surface. 

Proceeding along the lines in Ref. ^| we can now find 
the following approximate analytical solution for the den- 
sity of states spatially averaged over a unit cell C\ of the 
vortex lattice for arbitrary field directions: 



II. QUASIPARTICLE SPECTRUM IN THE 
VORTEX STATE 



Our calculation of the quasiparticle spectrum in 
the vortex state is based on quasiclassical Eilenberger 
theoryiiSiiii As we have shown in Ref. ^] there exists an 
accurate approximate method for the calculation of the 
density of states averaged over a unit cell of the vortex 
lattice (see section IV. A in Ref. 0) and we briefly repeat 
the essential equations here. We start from the Eilen- 
berger equation for the normal and anomalous compo- 
nent g and / of the quasiparticle propagator of a spin- 
singlet superconductoni^ii 6 . 
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where the momentum and frequency dependent function 
Pa takes the form 
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The w-function, also known as Dawson's integral, is re- 
lated to the complement of the Error function by 
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Here, A is the vector potential and is chosen to be A = 
— if x B, the vector v F — v F (k F ) denotes the Fermi 

velocity at the momentum point k F at the Fermi surface 
and e„ are the fermionic Matsubara frequencies. Eq. JJJ 
has to be supplemented by a normalization condition 15 
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+ f(r,k F ,ie n )f*(r,-k F ,ie n ) = l (2) 



For the spatial variation of the gap function A(r, k F ) = 
A(k F ) ■ tp\(x,y) we take the Abrikosov vortex lattice in 
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Furthermore z are normalized Matsubara frequencies of 
the following form 



V2e n 
\Vk F \ 



(9) 



where r}k F is proportional to the projection of the Fermi 
velocity into the complex plane perpendicular to the mag- 
netic field direction and given by 



K v f,i + iv F , 2 )\ 



(10) 
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Here, Vp,i and vf,2 are the components of the Fermi ve- 
locity in the plane perpendicular to the average magnetic 
field B. Using the flux quantization condition 



wilm u>2 
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where wilm L02 = \C'a\ is the size of a unit cell of the 
vortex lattice (if we assume u\ to be real) we can bring 
r]k F into the form 



r\k F = K v f,i + Wf,2) 



(12) 



If the magnetic field is not directed parallel to the c- 
axis of the uniaxial systems considered here, then the 
rotational symmetry around the magnetic field direction 
is broken and we have to take into account, that the 
isotropic vortex lattice no longer has to be the appro- 
priate ground state. If we use a variational ansatz for a 
distorted Abrikosov vortex lattice with 



Vk{x,v) = -0A(e T x,e T y) 



(13) 



we find that the Fermi velocity components in Eq. (|1(J|> 
have to be replaced by the scaled Fermi velocities^ 

vf,i — e T i>F,ij and vp,2 — s~ T vp.2 (14) 

The distortion parameter r has to be found by a mini- 
mization of the free energy at a given field strength and 
temperature. 

Within our approximation the free energy difference 
between superconducting and normal state is given by 
the expression 



n s ~n N = -(nT |A(fcF)|2 V^(^) 
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A derivation of this expression can be found in Appendix 
1X1 For each value of r the gap function in Eq. (|15|) has 
to be calculated from the gap equation 




|A(fc F )| s 



g{kF,ii n )V^zw{iz) — 1 
T 



-|A(fc F )|^ln- 



(16) 



FS 



sclf-consistently. We note that the variational param- 
eter r enters Eqs. (TL5j) and ((TBj) only via the quantity 
|%f( t )I- A determination of r requires to minimize 
Eq. i|15|) . where for each trial value of r Eq. i|16|) has 
to be solved self-consistently. Fortunately, we were able 



to find analytical results for r in some of the cases stud- 
ied below, such that this minimization procedure could 
be avoided in these cases. 

The field dependent total density of states averaged 
over the Fermi surface and a unit cell of the vortex lattice 
is found from Eq. |jfj| 



N(B, E) = ( Rc 
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We note that the direction of the Fermi velocity enters 
this equation only via the quantity |?7fe F | and that only 
the components perpendicular to the magnetic field play 
a role. This is where the Fermi surface structure comes 
into play. 



III. PARAMETRIZATION OF THE FERMI 
SURFACE 

In order to parametrize the different Fermi surfaces we 
have to find an appropriate parametrization of the Fermi 
wave vector kF- 

k F {d,v) = k x (-&,ip)e x + k y {8,(p)e y + k z (-&,<p)e z (18) 

Using this parametrization we can calculate the outward 
oriented normal unit vector n F 



\t v x t&\ 



(19) 



where t v and t$ are two orthogonal tangential vectors 
that are parallel to the coordinate lines of tp and i? 



tip = d v kF{$, f), itf = ip) 



(20) 



The Fermi velocity vector is always directed perpendicu- 
lar to the Fermi surface and thus can be written in general 
as 



(21) 



The Fermi surface of layered systems as for example 
the high T c -cuprates can be described as a distorted 
cylinder. The distortion is due to a small c-axis dis- 
persion. In this case we can parametrize the distorted 
cylindrical Fermi surface as 

k x (k c ,if) = (k a b + ^cos(cfc c )^ cosip, 
k y (k c , ip) = (k ab + ^ cos(cfc c )) sin ip, 



k z (k c ,ip) 



(22) 



with a dimcnsionless c-axis dispersion parameter e c . The 
parameter k c is running from —n/c to it /c where c is the 
lattice constant in c-axis direction. The polar angle ip 
varies in the interval [0, 2tt]. 
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Another Fermi surface that we will discuss is a half- 
torus. In magnesium diboride with its hexagonal crystal 
structure the 7r-band possesses a tubular Fermi surface 
structure that can be approximated by a half-torus living 
at the border of the hexagonal unit cell of the reciprocal 
latticed In this case we have the following parametriza- 
tion 

k x {d,(p) = kp ■ (y + cos •&) cosy, 
kyl'djip) — kp ■ (y + cosz?) siny, 
k z ($,(p) = fc F -sintf (23) 

where the parameter v denotes the ratio between the 
radius of the meridian and the radius of the longitudinal 
curve of the torus. Again tp can take values between and 
2ir while i? varies between ir/2 and 3tt/2 to parametrize 
the concave half-torus. The spherical Fermi surface can 
be viewed as a special case of the toroidal Fermi surface 
with parameter v — 0. 

If we assume a crystal with different effective masses 
in the c-axis direction and within the a6-plane 
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we find an elliptic Fermi surface with the following 
parametrization 

k x ($,tp) — fc F • cost? cosy, 

ky($,(p) — kp ■ cost? sin y, 

k z ($,tp) = k F ■ ^m c /m ab sin$ 



(25) 

Here, •& takes values between —tt/2 and ir/2. In Fig.^ we 
show the four different Fermi surfaces discussed above. 
From these parametrizations we can now calculate the 
Fermi velocity and its projection into the plane perpen- 
dicular to the magnetic field. We will consider the two 
cases of a magnetic field applied parallel to the c-axis di- 
rection of the crystal and the case of an in plane magnetic 
field. 

For the cylinder with c-axis dispersion e c we assume 
that e 2 c <C 1. In this case the variation of the in-plane 
components of the Fermi velocity due to the c-axis dis- 
persion can be neglected. Then we can write 



vp = vp(cos ip e x + sin y e y + e c sin ck c e z ) 



(26) 



For the half-torus and the sphere the Fermi velocity is 
found to be: 

vp = vp (cosy cos e x + siny cosi9 e y + sini9 e z ) (27) 

Writing vp^ — kp/m a b and vp, c = kp / \jm c m a b we find 
for the elliptical Fermi surface: 

v F = V k e k = (28) 
VF,ab(cos y cos "S e x + sin y cos •& e y ) + vp tC sin ■& e z 

The Fermi surface average Eq. (0 reads in these coor- 
dinates 
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for the cylindrical Fermi surface with f^<Cl and 
1 f 27T , , /" 37r/2 ,„i/ + cosi3 



fs 2tt J 



dd- 



tt/2 



TXV - 2 



(30) 



for the half-torus. The integrations for spherical and 
elliptical Fermi surface are obtained from Eq. I|3U|) for 
v = 0. 



A. Field in c-axis direction 

As pointed out above, for the calculation of the density 
of states in the vortex state using Eq. lfT7ll we only need 
to know the modulus of the quantity rjk F introduced in 
section [H] This quantity only contains the components 
of the Fermi velocity perpendicular to the external mag- 
netic field. For field directed in c-axis direction we there- 
fore only need the a6-plane components. As a result of 
the rotational symmetry of all the Fermi surfaces consid- 
ered we expect no angular dependence on the polar angle 
y. For the cylindrical Fermi surface we find no angular 
dependence at all 




(31) 



For the toroidal or the spherical Fermi surface a depen- 
dence on the azimuthal angle is found: 



\Vk F \ = fop |C0SI?| \ — 

V nc 

For the elliptic Fermi surface we have 



eB 

\Vk F | = flVF, ab |COS0| \J — 



(32) 



(33) 



B. Field in afo-plane direction 

If we consider a magnetic field that is not parallel to the 
c-axis direction, the rotational symmetry around the c- 
axis is broken and we find more complicated expressions 
for the Fermi velocity projection rjk F ■ In this case we 
can find the projection by a simple rotation, for example 
around the a-axis. If we carry out a rotation with angle 
7 we obtain a modified expression for the Fermi velocity 
in the rotated coordinate frame 

v F = v F<a e\ + (cos 7 v F , b - sin 7 vp^ c )e 2 

-I- (sin 7 vp tb + cos 7 v F ^)e 3 (34) 

Here the unit vector e*3 is chosen to point into the di- 
rection of the magnetic field and e\, e*2 span the plane 
perpendicular to the magnetic field. vp^ a , vpj, and vf jC 
are the components of the Fermi velocity along the a-, 
b- and c-axis of the crystal, respectively. If the magnetic 
field is applied along the a6-plane direction of the crystal, 
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the angle 7 equals 7r/2 and the Fermi velocity v' F in the 
rotated coordinate frame can be written as 



VF, a ei + (-VF,c)e 2 + VF,be3 



(35) 



As we have pointed out above, in this field direction the 
distortion of the vortex lattice has to be taken into ac- 
count which leads to a scaling of the components of the 
Fermi velocity (see Eq. (|14|)h In this way we find for the 
cylindrical Fermi surface 



\Vk F \ = hvF\ 1— \/e 2T cos 2 ip + e~ 2T e 2 sin ck c , (36) 
he 



for half-torus and sphere we have 
feB 



\Vk* 




h v F\l -i—\/e 2T cos 2 if cos 2 d + e~ 2r sin (37) 
he 



and for the elliptical Fermi surface we find 

feB I 

\Vk F I = W — y e2T v F ab cos 2 ip cos 2 ■& + e~ 2T v 2 F c sin 2 § 

(38) 

In this case it is useful to introduce a rescaling of the 
distortion parameter r via 



VF,ab 



VF.c 



(39) 



Then we can write 



\Vk F \ = fty/VF,abVF,c\J t"^ 21 cos 2 (p cos 2 + e~ 2T sin 1? 

(40) 

This equation has the same form as Eq. I|37|l with r being 
replaced by f and by y/VF,abVF,c- In the next section 
these expressions are used for evaluation of the Fermi 
surface integration and calculation of the spatially and 
momentum averaged density of states. 



IV. 



FERMI SURFACE INTEGRATION 



In order to calculate the quasiparticle spectrum, we 
have to integrate the spatially averaged density of states 
over the Fermi surface using Eq. I|17fl . In the following we 
are going to discuss the case of an isotropic s-wave gap 
and the two gap superconductor MgB2. Some results for 
a d-wave gap and a cylindrical Fermi surface can be found 
in our previous work^4 We note that for an isotropic 
s-wave gap the function Pa depends on momentum kp 
only via the quantity \r]k F \- In this case it is useful to 
introduce a characteristic Fermi surface function <7f(s) 
of the following form 



9f(s) 



N(0) J FS (2tt) 3 \hv F {k F )\ 



\Vk F 



(41) 
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FIG. 2: Averaged density of states in the vortex state for the 
cylindrical Fermi surface with magnetic field in c-axis direc- 
tion for three different magnetic fields B/B C 2 = 0.1, 0.3, and 
0.5. The inset shows the field dependence of the gap ampli- 
tude in units of Ao obtained numerically from Eq. 11611 (dots) 
and the approximation Eq. 148H (solid line). 



For convenience, we have introduced a parameter a con- 
taining the field dependence 



(42) 



For a given value of s, 3f(s) counts the number of states 
at the Fermi surface with a certain projection of the 
Fermi velocity into the plane perpendicular to the mag- 
netic field. Using the function <7_f(s) the density of states 
for a given value of a can be written as a single integra- 
tion over the variable s: 



N(B,E) 
with 



ds g_p(s) Re 



1 



v/l + P A (as,P) 



(43) 



P A (as,E) - 



4IAI 



r(l - \/Trz(as)w(iz(as))) (44) 



(a s) 2 

Here, z(as) is the normalized real axis frequency 

-i-v/2 



z(as) 



a s 



(E + 



(45) 



All information about the Fermi surface structure is con- 
tained in the function <7_f(s), which still depends on the 
direction of the magnetic field, but not on the magni- 
tude. Looking at Eq. I|43() #_f(s) can be interpreted as a 
weighting function. 



A. Field in c-axis direction 

The characteristic function gF{s) assumes a particu- 
larly simple form when the magnetic field is directed 
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FIG. 3: The characteristic functions gw[s) for the spherical 
and half-toroidal Fermi surface. The magnetic field is applied 
along the c-axis direction and we have chosen the ratio of the 
toroidal radii to be v = 4. 



along the crystal c-axis. For the cylindrical Fermi surface 
\r]k F | becomes angular independent and the characteristic 
function reduces to a <5-function: 

g F (s) = 5(s-l) (46) 

In this case the density of states is just given by 



JV cy i (B,E)= Re 



1 



y/1 + P A (a,E) 



(47) 



in agreement with Eq. (36) in Ref. 1J. This result is 
shown in Fig. [5] as a function of the quasiparticle energy 
E in units of the zero field gap value Ao = A(B = 0) 
for three different field strengths. For simplicity we have 
used a field dependence of the gap function 11 



A(S) = A(u/1- 



B 

B C 2 



(48) 



shown as the solid line in the inset of Fig. |2 This field 
dependence is a very good approximation to the field 
dependence obtained from the gap equation Eq. (|16l) nu- 
merically (solid squares in the inset) i24 The value of \rik F I 
at the upper critical field is found from the linearized gap 
equation Eq. (fT?)|l to be 



(49) 



Vk F \{B c2 ) = V2 7 A ~ 1.887A C 



with In 7 = 0.577215 being Euler's constant. 

As becomes apparent from Fig. |3 the peak-to-peak 
distance of the gap structure in the density of states 
increases with increasing field, while the gap itself de- 
creases. For this reason the gap structure seen in the 
density of states is not a good quantitative measure of 
the gap in the vortex state anymore. 

Apparently, the density of states in the general case 
Eq. i|43[) can be reduced to an integral over the density 



FIG. 4: Averaged density of states for the toroidal Fermi 
surface with magnetic field parallel to the c-axis direction 



of states for the cylindrical Fermi surface Eq. (|47|l : 



N(B,E) = dsg F (s) N cyi (as, E) (50) 
Jo 

This means that N(B, E) consists of a weighted average 
of the density of states for the cylindrical Fermi surface 
in Fig. |21 for different magnetic field strengths ranging 
from B = (corresponding to s = 0) up to a certain 
upper limit given by the maximum of \r]k F \ on the Fermi 
surface. 

For the toroidal Fermi surface the integral over the 
Fermi surface in Eq. (|41|l can be done analytically and 
for the characteristic function we find 



9f{s) = 



2vT 




for < s < 1 



else 



(51) 



The result for the spherical Fermi surface is obtained for 
v = 0. The same result is also found for the ellipti- 
cal Fermi surface, if vp in Eq. (|42|l is replaced by vp^. 
The weight functions for the half-torus and the spheri- 
cal Fermi surface are shown in Fig. |3J For the half-torus 
v = 4 (appropriate for MgB2) has been chosen. 3_f(s) for 
the half-torus starts with a finite value at s = 0. This 
has a topological reason: on the half-torus with field in 
c-axis direction there exist two lines of points at which 
the Fermi velocity is parallel to the magnetic field. At 
these points the projection of the Fermi velocity onto the 
plane perpendicular to the field vanishes. In contrast, 
for the spherical and elliptical Fermi surfaces there exist 
only two such points, the poles. Therefore gp(s) becomes 
at s — and increases linearly. For the cylindrical 
Fermi surface no such points exist and <?f(s) becomes 
for s ^ 1. 

In Fig.0]we show the density of states for the half-torus 
for different field strengths in c-axis direction. Because of 
the finite value of <?_f(s) for s = in this case there is large 
weight to the BCS singularity of the spectrum at B = 0. 
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FIG. 5: Averaged density of states for spherical and elliptical 
Fermi surface, respectively, with magnetic field parallel to the 
c-axis direction 



FIG. 6: The characteristic functions gp(s) for the cylindrical, 
spherical and half-toroidal Fermi surface. The magnetic field 
is applied in afe-plane direction. 



This results in the sharp flank at the gap edge seen in 
Fig. 0| even for higher magnetic fields. The position of 
the sharp flank closely follows the field dependence of the 
gap. The singularity itself is washed out by the averaging 
process. But still the most important contribution comes 
from the highest field, corresponding to s = 1 in Fig. [3J 
The characteristic function of the spherical Fermi sur- 
face increases linearly from s = 0. This results in a visi- 
ble break at the gap edge in the density of states but no 
longer in a significant flank, as shown in Fig. |SJ Again, 
the most important contribution comes from the highest 
fields (s = 1), leading to a peak-to-peak distance that in- 
creases with increasing magnetic field. The reduction of 
the BCS singularity is less pronounced than for the half- 
torus, because the weight near s — 1 in #f(s) (Fig. 0) is 
stronger for the sphere than for the half-torus. 

B. Field in a&-plane direction 

In order to calculate the characteristic function <7f(s) 
for field directed in a6-plane direction it is necessary to 
take into account the vortex lattice distortion via the 
parameter r. For each value of the magnetic field r has 
to be found by minimizing the free energy of the system 
as described in section [H] This makes the calculation of 
Sf(s) more involved in this case. However, even without 
knowing the value of r one can get already insight into the 
shape of the characteristic function <?f(s) by analyzing 
the extrema and saddle points of \rjk F \ in Eqs. (|36|l and 
l|37p. Via Eq. 1|41[1 these saddle points will lead to van 
Hove singularities in <7f(s). 

For the cylindrical Fermi surface we find from Eq. i|36|) 
that there is a minimum at |?7fc F |(<y9 = n/2,k c = 
0) = 0, a maximum at l^fc^Kv = 0, ck c — ir/2)/a — 
y / e 2T + e~ 2r el, and two saddle points at \r]k F \(^> — 
0, k c = 0)/a = e T and \r]k F \(<p = ir/2,ck c = ir/2)/a — 
e~ T e c leading to logarithmic singularities. In this case we 
can show analytically that r is minimized by e T = ^JT~ C 



for all field strengths. Therefore, the two saddle points 
collapse into one. The characteristic functions obtained 
after minimization of t are shown in Fig.Elfor field in ab- 
plane direction. The dotted line shows the result for the 
distorted cylinder. The logarithmic van Hove singularity 
due to the two saddle points is seen at s = yJT~ c ~ 0.4 and 
a step due to the maximum is found at s — \J2e c ~ 0.57. 

For the other three Fermi surfaces from Eq. (|37|1 we 
find a minimum at |?7/c F |(</9 = 7r / 2 , i? = 0) — and 
two extremal points at |?7fc F |($ = ir/2)/a = e~ T and 
\rjk F \((p = 0, i? = 0)/a = e T . The point at \rjk F \/a = e~ T 
for r > turns out to be an extended saddle point, lead- 
ing to a square-root singularity in gp(s). In the case of 
a spherical or elliptical Fermi surface r = or f = is 
found due to symmetry, respectively. As we have noted 
before, both the gap amplitude and the free energy func- 
tional depend on r only via \rjk F |- Also, as we have shown 
in Eq. I|40fl . we can bring |?7fe F | of the elliptic Fermi sur- 
face into the same form as \r]k F \ of the spherical Fermi 
surface with t substituted by f and vf substituted by 
y/VF,abVF,c- In this case the two extremal points again 
collapse into one and we just rediscover the same func- 
tion 9f(s) as for field in c-axis direction, as expected. 
For the half-torus the parameter r is nonzero and field 
dependent, however. Therefore, a square-root singular- 
ity is seen in Fig. |BJ at s = e~ T and a step due to the 
maximum at s — e T . Here, 3f(s) is shown for high fields 
close to B c2 where we find r = 0.375. 

In Figs. 0and|Hlwe show the corresponding densities 
of states for cylindrical and toroidal Fermi surface, re- 
spectively. For this field direction the densities of states 
for spherical (Fig. EJ, distorted cylindrical (Fig. , and 
toroidal (Fig. |HJ) Fermi surfaces appear to be much more 
similar to each other. This becomes clear from the simi- 
larity of the characteristic functions in Fig. H3 all curves 
increase linearly at low values of s. Again, the reason 
for this is topological: for field applied in a6-plane direc- 
tion in all three cases there exist a few singular points on 
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FIG. 7: Averaged density of states for the distorted cylinder 
with magnetic field in afc-plane direction 



LU 



2 
1.5 
1 

0.5 





B=0.1 Bfc 2 




ft - - - B=0.3 Bo2 






\ B=0.5 B c2 









-2 





E/A 



FIG. 8: Averaged density of states for the half-torus with 
magnetic field applied in afe-plane direction 



the Fermi surfaces at which the Fermi velocity is directed 
parallel to the magnetic field, in contrast to the case of 
field applied in c-axis direction. Accordingly, the results 
in Figs. [7| and |H1 resemble that of the spherical Fermi sur- 
face and we can identify again the characteristic break at 
the gap edge and the distinct peaks on each side of the 
gap that we have already recognized in Fig. [S] 



V. COMPARISON WITH A BROADENED BCS 
TYPE DENSITY OF STATES 

In analysis of experimental tunneling data on super- 
conductors in the vortex state a simple model of a BCS 
type density of states with a broadening parameter T of 
the following form has been used frequently^iiLi^ 



N(E) = Re 



E + iT 



y/(E + iT) 2 - A 2 



(52) 



for E > 0. Here, A and T have been used as field depen- 
dent fitting parameters. In this section we want to make 



FIG. 9: Spatially and momentum averaged density of states 
for B = 0.3 Bc2- The boxes show the data from the quasi- 
classical calculation (also seen in Fig. |2J and the solid line 
represents the broadened BCS spectrum that was fitted to 
these data using Eq. (1521 . 



a critical comparison of this model with our results and 
point out the limitations of this model. 

Looking at Figs. [3 0] and [SJit is immediately appar- 
ent that Eq. H52[) is certainly not able to reproduce these 
very different shapes of the curves, because the results in 
Figs. [21 El and [5] are shown for the same reduced mag- 
netic field and gap values (only the Fermi surface struc- 
ture was changed). Nevertheless one might ask whether 
these curves can be approximated by Eq. I|52[) with some 
effective broadening parameter T. In order to answer this 
question we have made least squares fits of Eq. 1|52[) to our 
results for the cylindrical Fermi surface Eq. (|47|l shown 
in Fig.0 In Fig. ©we show the best fit for B = 0.3B c2 . 
It appears that the width of the peaks at the gap edge 
and the high value of the density of states at small ener- 
gies cannot be described simultaneously very well. The 
main reason for this failure is due to the presence of a 
high amount of vortex core states that dominate at low 
energies and cannot be captured very well by Eq. (|52|l . 

In Fig. we show the values of the two parameters 
A and T as a function of magnetic field found from our 
fits (solid squares). In the upper panel we compare A(£?) 
with the actual field dependence of the gap shown as the 
solid line. This comparison shows that the gap value ex- 
tracted from Eq. (|52|l does not very accurately reproduce 
the actual gap value and care should be taken when in- 
terpreting gap values extracted this way. The lower panel 
in Fig. I1UI shows the field dependence of the broadening 
parameter T. It is clear that this quantity increases with 
increasing magnetic field. For illustration we are also 
showing the 'average Doppler shift' \rjk F |/v2 as the solid 
line. This quantity represents the average energy broad- 
ening due to the local Doppler shift of the energies of the 
excited quasiparticles in the presence of the supercur- 
rents around the vortices. The comparison shows that T 
roughly measures this average Doppler shift, particularly 
at higher magnetic fields. 
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FIG. 10: The upper panel shows the gap amplitude that was 
extracted from the fitting process (boxes) in comparison with 
the field dependence of the gap function that was used to 
calculate the quasiclassical spectra (solid line). The lower 
panel shows the variation of the fitting parameter F (boxes) 
as a function of the magnitude of the magnetic field compared 
with the average Doppler shift \rjk F \/\/2 (solid lines). 



Instead of using Eq. I|52l) for an extraction of the mag- 
netic field dependence of the gap from experimental data 
it would be better to use the theory presented here. This 
requires a priori knowledge of the Fermi surface struc- 
ture, however. For the special but important case of the 
cylindrical Fermi surface with field in c-axis direction the 
calculation simplifies and we can suggest an improved 
version of Eq. (|52|l : first we note that Eq. I|47(l depends 
on magnetic field only through the field dependence of the 
gap function A(B) and the quantity a oc \fB. Therefore, 
Eq. 147|) can be brought into the following convenient 
form: 



N(B, E) = Re 



1 



y/l + P A (B,E) 



(53) 



with 



P A {B,E) = 2^ (^1 ) (l + iV^xw{x + i0+)) (54) 



B V A 



and 



B E 



x = iz — \ — 

BA 



(55) 



where we have introduced a characteristic magnetic field 



B 



2c 



hvi 



A 2 



(56) 



such that \ % = 75rT ■ In Eq. I|53l) all material depen- 



dent quantities including the Fermi surface structure ap- 
pear lumped into the single parameter B, which makes 
Eq. (|53|l particularly useful for fitting of experimental 
data. Such a fitting could proceed as follows: first, the 
zero field value Ao has to be extracted from a usual fit 
to zero field experimental data. Then, Eq. (|53|l can be 
fitted to finite field data using B as a field independent 
and A(S)/Ao as a field dependent fitting parameter. For 
this fitting procedure a useful approximation of the ex- 
pression in Eq. (|54|l containing Dawson's integral for real 
values of x is given by 



1 + i\firxw(x + i0 + ) 



1 



0.2x 4 



1 



0Ax e 



(57) 



This approximation is better than 2% and reproduces the 
behavior of the left hand side in the limit x — > to order 
x and in the limit x — * oo to order 1/x 2 . 

In order to get some feeling for the parameter B, we 
can relate it to the upper critical field B c2 using the lin- 
earized gap equation. At B C 2 we findi^ 



a(B c2 ) = 
and therefore 
B 

B C 2 



27A0 exp 



In 



\Vk f 



FS 



(58) 



exp {2 (in 1^) J 
a(B C 2) ) 7 



(59) 



For the cylindrical Fermi surface with field in c-axis di- 
rection considered here ^ VkF ■ = 1 and therefore we have 
B = B c2 h = 0.561B c2 . 

VI. APPLICATION TO MGB 2 

In this section we want to apply our findings to the sit- 
uation in MgB 2 . It has been established recently that this 
compound belongs to the rare case of a superconductor 
possessing two different sized gaps on two different parts 
of the Fermi surface: a large gap living on the cylin- 
drical er-bands and a small gap on the three-dimensional 
Tr-bandsiiSiSSiSii 2 ^ The structure of the 7r-band Fermi sur- 
face can be approximated by a half-torus as pointed out 
111 Ref. E3 The quasiparticle excitations in the vortex 
state of MgB2 have been studied recently by tunneling 



spectroscop; 



.4.5.6.22 



and the field dependence of the su- 



perconducting gap was extractedi 4 ' 5 ' 6 

The normalized total density of states of a two band 
superconductor is given by a weighted average of the two 
partial densities of states for each of the two bands: 

N(B, £) = (!- w„)N a (B, E) + w n N„(B, E) (60) 



Here, is the weight of the 7r-band density of states and 
has been calculated by band structure calculation: 



19,23 
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FIG. 11: The two gaps of MgB2 as a function of magnetic 
field. The solid squares show the maxima of the peaks that 
result from the a band density of states, the triangles mark 
the flanks of the smaller gaps, also extracted from the quasi- 
particle spectrum. 



be w-n = 0.577. For the partial densities of states we can 
take our results from section llVl for the distorted cylinder 
with the distortion parameter e c = 0.163 and the half- 
torus v — 4. The two different gap amplitudes have to 
be found from a solution of the two by two gap equation 
in the vortex state: 

A Q (r) = 7rT^ ]T \ a ' a ' (f a '(r,k Fa ,ie n )) Fs (61) 

a' |e„|<aj c 

where we assumed an isotropic pairing interaction X a,a 
with a the band index: a = a, n. 

In Fig. II II the two gap amplitudes are shown as func- 
tion of the applied magnetic field for the set of parameters 
used in Ref. [l3- The resulting density of states for field 
in c-axis direction is shown in Fig. 1121 From this figure 
it becomes apparent that the peak at the gap edge of the 
small gap is much more rapidly suppressed as the field is 
increased. This is due to the topology of the 7r-band: as 
we have shown in the previous section its half-torus shape 
leads to a much more rapid suppression of the peak at 
the gap edge than for the cylindrical Fermi surface of the 
(T-band. Experimentally, this effect has been noted by 
Gonnelli et ali^ Already at comparatively low fields of 
about 1 Tesla no apparent structure of the 7r-band gap 
could be observed anymore, while the structure of the a- 
band gap remains visible to much higher fields of about 
4-5 Tesla. 

For comparison in Fig. 1111 we are showing the position 
of the peak in the cr-band density of states as a function 
of magnetic field as the solid squares. The position of the 
sharp flanks in the 7r-band density of states is shown as 
the triangles. Clearly, the structures in the densities of 
states are behaving quite differently in the two bands due 
to their different topologies. This effect should be taken 
into account when analyzing the experimental data. 
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FIG. 12: Averaged density of states for MgE>2 with magnetic 
field parallel to the c-axis direction of the crystal lattice, cal- 
culated as a weighted sum of the averaged density of states of 
a toroidal Fermi surface and a cylindrical Fermi surface using 
Eq. I6UI . Results are shown for different magnetic fields. 



VII. CONCLUSIONS 

We studied the influence of Fermi surface topology on 
the density of states in the vortex state of type II super- 
conductors for the four Fermi surface structures shown 
in Fig. ^ The topology of the Fermi surface takes influ- 
ence on the density of states, because the direction of the 
Fermi velocity with respect to the local supercurrent flow 
around the vortices leads to a change of the excitation 
energy of the quasiparticles. We saw that the density of 
states behaves quite differently for cylindrical, spherical 
and toroidal Fermi surfaces. The field dependence and 
shape of the curves shows characteristic features related 
to the topology of the Fermi surface in question. We 
showed that these features can be understood in terms 
of characteristic Fermi surface functions. A particularly 
important role is played by the number of points on the 
Fermi surface at which the Fermi velocity is directed par- 
allel to the external magnetic field. 

We compared our results with the simple model of a 
broadened BCS type density of states, that has been used 
frequently in the past for analysis of experimental data. 
This comparison showed that the contribution coming 
from vortex core states is underestimated in this simple 
model. This leads to inaccuracies in the gap values ex- 
tracted from the broadened BCS model. For the special 
case of a cylindrical Fermi surface with field along the c- 
axis direction, we suggested an improved formula for the 
density of states in the vortex state, in which all material 
properties are lumped into a single parameter. This new 
formula does not possess the limitations of the broadened 
BCS model and can be used for fits to experimental data. 

We applied our results to the case of the two gap su- 
perconductor MgB2. This case is particularly interesting, 
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because the two Fermi surfaces related to the two gaps 
possess completely different topology. We demonstrated 
that this leads to very different field dependencies of the 
partial densities of states, resembling recent observations 
on this compound. 
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APPENDIX A: CALCULATION OF THE FREE 
ENERGY 

In this Appendix we provide a derivation of Eq. I|15|) 
for the free energy difference within the method we are 
using here. The difference of the free energy between the 
superconducting state and the normal state for the case 
of an even parity superconductor can be calculated by 
the coupling constant integration method^ 



— = I — {XH int )\ 



(Al) 



where Hi„ t denotes the interaction Hamiltonian. For a 
multiband superconductor this equation can be expressed 
in terms of the gap function A" (r , kp ; A) and the anoma- 
lous Eilenberger propagator f a (r, hp, ie n ; A) correspond- 
ing to a reduced pairing interaction A • V aa {kp, k' F ) 

- 1 dX r ^ 



t 



xnT f a (r,k F ,ie n ;X)\ (A2) 

|e„|<w c / 



Here, a is the band index and A Q (r, kp; A) is solution of 
the renormalized gap equation 



A Q (r,fc F ;A) 



(A3) 



X J2/v aa '(kp,k'p)TTT f a '(r,k'p,ie n ;X)\ 

a' \ |e„|<w c / FS' a , 

Here, /" (r,k' F ,ie n ',X) is a solution of the Eilenberger 
equation Eq. in the presence of the A-dependent 
gap function A Q (r, kp;X) and depends on the param- 
eter A only implicitly via A". The integration over A 
in Eq. (|A2|) can be substituted by an integration over 
the pairing potential, a usual technique as described in 
textbooks^ Generalizing this to the case of an inhomo- 
geneous superconductor we introduce a function x(X) via 



Then, f a can be viewed as a function of x, calculated for 
a reduced gap iA°, and the integration over A can be 
substituted by an integration over x. In order to do so, 
we insert Eq. (|A4jl into Eqs. (|X2|) and (JX3J). Taking the 
derivative of Eq. (|A3|) with respect to x we find 



(A5) 



t dX ^ 

- — A«(f,k F ) = A Q (r,M- 
A dx 



Xj2\V aa '(kF,k'p)7rT ^-( r ^' ie « 



|e»|<u 



FS' , 



Inserting the hermitian conjugate of Eq. I|A5|I into 
Eq. I)A2|I and using the fact that the pairing interaction 
is hermitian. i.e. 

(A6) 
(A7) 

A"ir./, : F )j' f a (r, k F ,i€ n ;x) 

x——{r, kF,ie n ;x) 
ax 




FS a 



A partial integration of the last term with respect to x 
finally yields 



N 



C A 



d 2 r 

|Ca| 



|e„|<^c 



t 



f a (r,k F ,ie n )-2 / dxf a {f,k F ,ie n ;x))j) (A8) 

This equation still holds generally for inhomogeneous 
multiband superconductors. 

We can now employ the approximation procedure in- 
troduced in Ref. [lj and evaluate the spatial average in 
Eq. I|A8(1 . For a single band superconductor we find 

t -. - 



1 



/ d 2 r 


A(f,k F ) 


lc A 





\Ca\ 

|A(M| 2 



^TTZ w(iz) 



1 + x 2 PA(k F , ie r , 



(A9) 



a; (A) • A Q (f, k F ) = A°(r, k F ; A) 



(A4) 



The first term in Eq. i|A8(l is obtained for x = 1, of course. 
In this expression z denotes the normalized Matsubara 
frequencies as introduced in Eq. Q. Using the same 
method Eq. (|16|l can be derived from Eq. l)A3Jl for A = 1 
averaging over a unit cell of the vortex lattice and elimi- 
nating the pairing interaction V using the linearized gap 
equation at T = T c . Finally, the integration over x in 
Eq. I|A8(1 can be carried out analytically using Eq. (|A9|) 
and we obtain the result in Eq. (|15|l . 
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